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Abstract. A parametric study of Drude and Lorentz models performances in 
maximizing near-field radiative heat transfer between two semi-infinite planes 
separated by nanometric distances at room temperature is presented in this paper. 
Optimal parameters of these models that provide optical properties maximizing the 
radiative heat flux are reported and compared to real materials usually considered 
in similar studies, silicon carbide and heavily doped silicon in this case. Results 
are obtained by exact and approximate (in the extreme near-field regime and the 
electrostatic limit hypothesis) calculations. The two methods are compared in terms 
of accuracy and CPU resources consumption. Their differences are explained according 
to a mesoscopic description of near-field radiative heat transfer. Finally, the frequently 
assumed hypothesis which states a maximal radiative heat transfer when the two 
semi-infinite planes are of identical materials is numerically confirmed. Its subsequent 
practical constraints are then discussed. 
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It has been shown in the late 1960 ^ [2] that the radiative heat flux (RHF) exchanged 
by two media in the near-field (NF), i.e. when these media are separated by very small 
distances (smaller than the thermal radiation characteristic wavelength At = ^) could 
exceed by several orders of magnitude the black body limit. This topic has then received 
an increasing attention until its recent experimental verifications [21 SJ . Experiments 
exclusively focused on asymmetric configurations such as plane-tip or plane-sphere 
configurations. On the other hand, the symmetric plane-plane configuration, potentially 
useful for various applications such as the cooling of high fiux density electronic devices 
[6] or thermo-photovoltaic (TPV) conversion of radiative energy [7], has been thoroughly 
investigated from a theoretical point of view by several groups ^ These 
theoretical works mainly addressed dielectrics, usually silicon carbide (SiC) [HI [9] 
which surface phonon-polaritons highly contribute to the NF RHF increase. They 
also considered materials which support plasmon-polaritons in the wavelength range 
of thermal radiation at room temperature such as tungsten [121 E] or heavily doped 
silicon (HD-Si) [iQlttS]. 

In the present numerical work, hypothetical materials modeled by local Drude and 
Lorentz models are considered. The aim of this work is to find the sets of parameters 
of these models that possibly maximize the RHF between two semi-infinite planes 
of identical materials separated by a nonometric gap at room temperature. For 
this purpose, we calculate the exchanged RHF between the two media while varying 
the different parameters in a wide range. Exact and approximate calculations are 
performed and their accuracy /resources consumption ratio compared. Then, the 
optimal hypothetical material performances are compared to those of usually considered 
materials, SiC and HD-Si for instance. Finally, the infiuence of small discrepancies 
between the optical properties of the two planes on the exchanged RHF is discussed. 

2. Formalism 

The two methods used to calculate NF RHF between two semi-infinite planes and to 
obtain results presented later in this paper are briefiy reminded in this section. 

2.1. Exact calculation 

Consider two semi-infinite planes 1 and 2 separated by a gap of thickness 6 (Figure 
[T]) and characterized by their dielectric functions and temperatures (ei,Ti) and (e2,T2) 
respectively. The total RHF density exchanged by the two media is given by [T4] : 

Q Qprop H" Qevan (l) 
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Figure 1: Two semi-infinite planes separated by a distance 5. 
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are the contributions of propagative and evanescent waves respectively. 
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is the mean energy of a Planck oscillator at a temperature T. r\ - are Fresnel refiection 
coefficients for an i-polarized wave {i G {s^p}) propagating from medium 3 to medium 
j. 73 is the wave vector normal component in medium 3 and is given by: 




73 = y(^-; (5) 

where q is the component of the wave vector parallel to the interfaces. Expressions [2] 
and |3] express the total heat ffux as the sum of the energy of different existing oscillators 
at a temperature T, transported by different modes (cL;,q). It is worth noting that for 
Q > ^5 73 is imaginary, the corresponding waves are evanescent and their magnitude 
decreases when going away from the surface. Corresponding modes are surface waves 
modes. 

To obtain the total heat ffux, a double integration over all modes (cj, q) is to be made. 
Its calculation may prove to be very resource-consuming since the cutoff wave vector qc 
for the integral over q is not known a priori. Different authors have proposed different 
approximations for the cutoff wave vector : qc = 1/5 |T3l [I5], qc ^^1^^ + {^/^Y [ISj 
and qc = 1/a [17j where a denotes the lattice constant of the considered material. In 
this work, = 50/5 is adopted. 



2.2. Approximate calculation 

Recently, Rousseau et al. [T3l[l8l|9] derived, under few simple conditions, an asymptotic 
expression of the NF RHF p-polarized evanescent contribution. This contribution is 
considered for two reasons. First, it dominates the other contributions in extreme near- 
field regime for dielectrics and some other materials such as HD-Si for instance. Second, 
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its exact calculation is the most resource-consuming due to the unknown and eventually 
large cutoff wave vector in some situations. 

First, they started, when considering a small temperature difference ST between the 
two planes, by defining a radiative NF exchange coefficient h : 

m T) 



h 



5T 



which can be written as the sum of two coefficients hprop and hevan 



(6) 

corresponding to the 

propagative and evanescent contributions respectively. Let's focus on the z-polarized 
{i G {p, s}) monochromatic evanescent contribution to the radiative transfer exchange 
coefficient which is given by : 
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where Tl^ani'^zj^ ^) (i^, q) mode transmission probabihty from medium 1 to medium 2 
dSmg and 



1 hu huj^ 



1 



2 sinh(^), 



(9) 
(10) 



TkTAn^c^ 

is proportional to the Planck function derivative. If we consider the electrostatic regime, 
i.e. Q Qo = p-polarization Fresnel coefficients become independent of q since they 
tend toward = Then, we can show [S] that /letiom prevailing in our case, may be 
written as : 
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evaluation), h^{u) = j-^z^^ = 
NF heat flux is then given by : 
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Therefore, the heat flux calculation is reduced to a simple integral evaluation and the 
problem of the cutoff wave vector is apparently resolved. Given the assumed hypotheses 
in order to obtain expressions [TT] and [12} a verification with an exact calculation of 
results obtained by this method might be necessary. 



3. Optimization State of the art 



Different groups have already tackled the question of maximizing the NF radiative heat 
transfer, for plane-plane configuration in particular. Zhuoming Zhang's group of Georgia 
Tech. has been particularly prolific. First, Basu et al. [21] led a theoretical parametric 
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study of radiative transfer between two semi-infinite planes of HD-Si. This material was 
considered because of its interesting optical properties that can be controlled through 
the doping level [211 ESI [23]. In fact, its dielectric permittivity is modeled by a Drude 
model where the doping concentration controls both of the plasma frequency uUp and 
the damping coefficient F. They observed that the RHF spectrum presents a peak 
around the plasma frequency and a blue-shift of the peak position when the doping 
concentration increases. They also noted that the total exchanged RHF increases with 
doping until a maximum that depends on temperature and ujp. At room temperature, 
this optimum is observed for a doping concentration between 10^^ and 10^° (cm~^). Let 
us note that similar results, obtained by a different approach, have been reported for 
HD-Si by \13\. Finally, they considered two planes with different doping concentrations 
and tend towards the conclusion that the maximal RHF is obtained for identical media. 
Then, in another work [24j, they went beyond the particular case of HD-Si by considering 
two identical semi-infinite planes of a completely fictive material. They found that the 
dielectric permittivity maximizing the exchanged RHF can be written e = —1 + iS with 
Im{e) = (5 <C 1. It is worth noting here that this form of e underlies a hypothesis 
of a non-dispersive medium. At the same time, Wang et al. [17j considered less 
restrictive situations and generalized first results previously obtained for HD-Si to other 
real materials (SiC, MgO) and fictive materials modeled by Drude and Lorentz models. 
For Drude model, control parameters are uj^ and F and the high frequency limit of 
the dielectric permittivity e^o- Lorentz model has an additional parameter cjq which 
corresponds to the frequency of transverse optical phonons. Authors make the following 
general conclusions : (1) Drude model leads to higher values of maximal RHF than 
Lorentz model. For this reason, Lorentz model presents its highest performances when 
cjo = 0, i.e. when it is equivalent to Drude model. That is why we focus on Drude 
model in the following points. (2) For Drude model : (2-1) Lower values of lead to 
the highest values of maximal RHF. These values are the closest to the condition given 
by [24j and previously presented. (2-2) At room temperature, a maximum of RHF is 
observed for uj^ 10^"^ (rad.s"-^) and V /uj^ 0, L The position of this maximum is 
strongly T-dependent. In addition, the maximum is realized by a compromise between 
the peak width (controlled by F) and the peak position (controlled by Up). 
More recently, several authors exploited graphene features to enhance NF RHF. 
Graphene presents palsmon-polaritons in the terahertz domain which makes it 
particularly interesting for radiative heat transfer around room temperature [25]. 
Besides, more than HD-Si, its optical properties can be tuned with doping level or 
chemical potential. Finally, graphene dielectric function is non-local, i.e. its dielectric 
permittivity in general, and its plasma frequency in particular, depend on the wave 
vector. Therefore, it presents a big variety of resonant modes which may allow to 
consider their coupling with other materials resonant modes. These authors showed 
that a thin film of graphene deposited on a dielectric that does not support surface 
phonon-polaritons leads to an enhancement of the exchanged NF RHF between two 
semi-infinite planes of the same graphene-covered material by three and almost four 
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orders of magnitude. However, this enhancement decreases rapidly with temperature 
and is spectacular only for temperatures lower than room temperature. Another group 
from rinstitut d'Optique of Paris [26], showed for a plane-plane system of SiC, that a 
thin film of graphene on the surface of one of the two planes leads to additional peaks 
in the spectrum of the local density of states due to the coupling of graphene modes 
with those of SiC. These modes contribute to the increase of exchanged NF RHF. 
Shall we here emphasize practical potential of graphene as a selective emitter for NF 
TPV devices. Indeed, the possibility to tune graphene plasmon-polariton resonance 
frequencies would allow their adjustment to the band gap of different photovoltaic 
converters. Messina et al. actually demonstrated [27j for a TPV device composed of a 
boron nitride emitter (at Tg = 450 K) and an indium antimonide cell, that a graphene 
film with a chemical potential of 0.5 eV on the surface of the cell leads to an increase of 
the maximal efficiency of the system by a factor two to reach r] 20% and an increase of 
output power by almost one order of magnitude. Higher performances corresponding to 
higher operating temperatures in the range [600, 1200] K have been recently presented 
by another group of the MIT [28j who considered a slightly different system where 
graphene plays the role of a selective emitter. Prior works had already considered NF 
TPV devices based on metallic selective emitters such as tungsten [I2[|8l|7] but graphene 
seems to monopolize the community recent attention due to the diversity of potential 
applications it makes possible thanks to the "ffexibility" of its surface modes and optical 
properties. 

4. Results 

Formalisms presented in the first section are used to calculate the exchanged NF RHF 
between two semi-infinite planes separated by a distance 6 = 10 nm. Planes dielectric 
functions are modeled by local Drude and Lorentz models usually adopted to describe 
real materials (gap thickness considered here is much larger than non-local phenomena 
onset distances [29l[30]). Calculations are made for both identical and different planes 
around 300 K while varying models parameters in their usual variation ranges with three 
main goals in mind : (1) For identical planes : to determine optical properties, fictive 
in this case, that would maximize NF RHF in order to guide, for a given application, 
the choice of a real material to use or the design of meta-materials (2) For different 
planes : to verify the hypothesis which states that the maximal RHF is obtained when 
the two planes materials are identical (3) To compare the accuracy and the resource- 
consumption cost of the exact and approximate methods. 
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4.1. Drude model 

We remind the expression of the dielectric permittivity in this model : 

e{u) = 600 - (13) 
where e^o is the high frequency hmit of the dielectric permittivity, ujp the plasma 
frequency and F the damping coefficient. 



4.1.1. Identical media 

First, we calculate NF RHF between two identical semi-infinite planes modeled by 
Drude model with = 1. (^p and F are varied within the ranges [10^^, 10^^] and 
[10~^ X cjp, 10 X ujp] respectively which cover these parameters ranges for HD-Si. Some 
values of these parameters for HD-Si with doping concentration around 10^^ cm~^ are 
given in Table [T] Media 1 and 2 are considered at 300 K and 299 K respectively. 



N° 


Doping type 


Concentration xlO ^^{cm ^) 


^00 


LOp X IQ-^"^ (rad.s-1) 


v_ 


1 


Si:B 


27 


11.8 


16 


6.7 X 10-1 


2 


Si:B 


6.7 


11.8 


8.3 


1.7 X 10-1 


3 


Si:P 


10 


11.8 


9.7 


5.1 X 10-1 


4 


Si:P 


5.3 


11.8 


7.28 


10-1 


5 


Si:P 


1.6 


11.8 


4 


1.3 X 10-1 


6 


Si:P 


0.52 


11.8 


2.3 


6 X 10-2 



Table 1: Drude model parameters for the dielectric permittivity of p and n-type HD-Si 
with bore (Si:B) and phosphorus (Si:P) respectively [31] . 



Plasma frequency and damping effects 

We present in figure [2] the normalized RHF exchanged between the two media. Plotted 
results are obtained by both exact calculation (Figure 2a) and asymptotic calculation in 
the case of extreme NF with the electrostatic limit approximation (Figure 2b). Only the 
dominating p-polarization is presented here. First, we can note the actual existence of a 
maximum (See Table [2] for its value and coordinates.). Beyond the maximum position, 
these figures reveal the RHF sensitivity to the different parameters. In fact, we observe 
that a relative variation between 22% and 27% for Up and of about 50% for F gives 
values of the fiux larger than 0.95 x qmax- These parameters values admissible variations 
to keep high fiux values are slightly larger than those reported in literature[10j. 



High frequency limit of the dielectric permittivity effect 
Similar calculation results are presented in Figure [3] for e 
and 600 = 10 (figures 3d and 3c). 
maximal fiux value when e^o increases. In fact, lower values of e^o lead to lower values of 



= 5 (Figures 3b and 3a) 
We observe as Basu et al.[10j, a decrease of the 
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Figure 2: Normalized NF RHF between two identical semi-infinite planes. The dielectric 
permittivity is given by Drude model with = 1. Results are obtained by both exact 
calculation (a) and asymptotic one (b). 



Method 


^00 




V_ 


qmax (W.m-2) 






t (s) 


E 
A 


1 
1 


1.51 
1.05 


0.17 
0.24 


229336 
229208 


500 
1000 


100 
100 


32921 
23 


E 
A 


5 
5 


2.51 
0.79 


3.7 X 10-2 
0.11 


78656 
78676 


500 
1000 


100 
100 


14864 
13 


E 
A 


10 
10 


3.47 
0.76 


1.51 X 10-2 
6.91 X 10-2 


42123 
43128 


400 
1000 


100 
100 


13952 
13 


E 
A 


20 
20 


4.57 
0.72 


2.29 X 10-3 
3.71 X 10-2 


24269 
22621 


400 
1000 


100 
100 


16961 
13 



Table 2: Drude model parameters maximizing the exchanged NF RHF between two 
semi-infinite planes separated by a gap of thickness 5 = 10 nm for different values of 
600- Values are obtained by both exact (E) and approximate (A) calculations, n^^^ and 
'^ujp.v are mesh points numbers for {uj^q) modes and for control parameters {ujp and F) 
respectively, t is CPU calculation time to obtain the corresponding figures. 



Re{e) = — R^{ ^2^-^^Yuj ) which are the closet to fit Basu et al. condition to maximize 
the NF RHFL24J, i^e. 'i?e(e) = -1. 



Exact versus approximate calculation 

Maximal values of the RHF obtained by both methods are almost the same with a 
relative error around 10~^ (see Table [2]). In the case 600 = 1, maxima are realized for 
J) = (1.51 X 10^^ rarf.5-^ 1.7 x IQ-^) and (cj^, J) = (1.05 x lO^^ rad.s-^, 0.24 x 
10~^) with exact and approximate calculations respectively. The relative error on 
positions is quite important, up to 27% and 82% for uj^ and F respectively. An exact 
calculation of the fiux value corresponding to approximate optimal parameters is 30% 
lower than the actual maximal flux value. This discrepancy on optimal parameters given 
by both methods increases with eoo- Let us note however the resource-consumption gain 
made by the use of the asymptotic approximation : figure [2a| (exact) was obtained in 
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Figure 3: Normalized NF RHF between two semi-infinite planes modeled by Drude 
model for = 5 : (a - b) and 6oo = 10 : (c - d) obtained by exact (left column) and 
approximate (right column) calculations. 



32921 (s) versus 23 (s) for figure 2b (approximate), i.e. a ratio of almost 1500 between 
the two. This ratio particularly depends on the parallel wave vector mesh resolution 
and increases rapidly with it. Calculations were made on an Intel® Xeon® E5620 @ 
2.40GHz, with 12288 Kb of cache and 4 Go of RAM memory. 
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Figure 4: Transmission probability of p-polarized evanescent modes r^yani^^ ^) defined 
in equation |8] for Drude model with = 1. Two cases are considered : (a) the exact 
optimum given in (Table [2| line 1) and (b) the approximate one given in (Table [2| line 

2). 



In order to understand the origin of the discrepancy between the two methods, we plot 
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in figure [4] each p-polarized evanescent mode (cc;,q) transmission coeflficient T^yani^^^) 



as defined in equation^ We consider the exact optimum (Table ^ fine 1, Figure 4a) 



and the approximate one (Table [2| line 2, Figure 4b). Several observations can be 
made : (1) The approximate optimum presents a high transmission coefficient (red and 
yellow areas) for more numerous modes than the exact one. (2) This higher number 
of transmitted modes is more pronounced for modes with a large wave vector parallel 
component q. (3) The exact optimum presents less transmitted modes but at a higher 
circular frequency. 

According to these observations, it is obvious that the discrepancy between the optimum 
position given by both methods is due to the fact that the electrostatic approximation 
ignores modes with low q. The optimum position shift in the approximate approach 
also induces a decrease in each mode mean energy. This decrease is compensated in 
the overall ffux density by a larger modes number. In order to accurately estimate each 
mode contribution, the value of the integrand (q/qo) x r^^ani^^ Q) ^^^^ Q which 

appears in equation [8] is more relevant than the mere modes transmission probability. 
It is plotted in figure [5| 
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Figure 5: The integrand (q/qo) x r^^ani^^Q) ^^^^ ^he wave vector parallel 

component in equation [8] for Drude model with 6oo = 1- The same cases as in figure [4] are 
considered : (a) the exact optimum given in (table [2| line 1) and (b) the approximate 
one given in (table [2| line 2). 



It appears through this figure that the weight of high wave number modes is dominating. 
These modes, for both exact and approximate optima, lay in the same q/qo range, 
q/qo G [200, 600] in this case, but for slightly different circular frequencies however which 
may explain the small relative error on ffux density values obtained by both methods. 
Finally, ffgure [5] allows an accurate calculus of the cutoff wave vector value. If we deffne qc 
as the largest wave number verifying (qc/qo) x rP^ani^p, qc) = \ [{q/qo) x rl^aA^P^ ^)\max^ 
we obtain (qc/qo) = 597 and (qc/qo) = 768 which leads to qc — 3/S and qc — 2.7/6 for 
the exact and the approximate calculation respectively. The cutoff wave vector is hence 
of the order of 1/5 and was actually overestimated in our first calculations. 
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Case of heavily doped silicon at 300 K 

Figure [6] presents similar results for e^o = e^o^si = 11-8. The aim here is to determine 
whether HD-Si, previously considered by several authors [lOl [T3] to maximize NF RHF, 
is well adapted to this task. For this reason, parameters values corresponding to HD- 
Si are represented on the same figure by crosses. Previously reported results [101 [I3] 
stating a maximal fiux for a doping concentration between 10^^ and 10^^ (cm~^) are 
more likely to be confirmed. Besides, we can state according to this figure that HD-Si 
around 10^^ (cm~^) is a good candidate to NF RHF maximization at room temperature 
since it allows to reach almost 0.9 x q^y^ax that can be obtained with a Drude model with 
6oo = 11.8 (we obviously assume that e^o is a parameter that can hardly be varied). 
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Figure 6: Normalized NF RHF for Cqo = Goo,si = 11-8 obtained by exact calculations. 
Crosses represent p and n-type HD-Si at different doping concentrations (see Table [l| 
for parameters values of the different points). 



4.1.2. Non-identical media 

The only change considered in this paragraph lies in the fact that exchanging semi- 
infinite planes dielectric functions are not identical while they are still modeled by 
a Drude model. We are aiming to a double objective : (1) verify the statement of 
maximal fiux for identical media due to a more efiicient coupling of identical modes 
supported by the same materials (2) See in what extent, a more or less important 
diflFerence in optical properties of considered materials affects the exchanged NF RHF. 
This second objective has an obvious applied interest since real materials eventually used 
in a particular application are never exactly identical. For this sake, we consider the two 
media around 300 K separated by 5 = 10 nm. We also consider, without generality loss, 
6oo = 1. Medium 1 parameters are fixed to optimal values previously obtained (Table 
[2| line 1). Control parameters are then the second medium Drude model parameters, 
i.e. (jjp^2 and Figure [7| presents the normalized exchanged NF RHF as a function of 
plasma frequencies ratio ujp^2/^p,i and damping factors ratio r2/ri. 
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Figure 7: Normalized NF RHF exchanged by two semi-infinite planes of different optical 
properties modeled by Drude model as function of plasma frequencies ratio and damping 
factors ratio obtained by exact calculation. The figure is obtained by exact calculations. 

The maximum is actually realized for {(jJp^2/^p,i^^2/^i) = (l^l), i-^. for identical 
media. Besides, the fiux value is more sensitive to Up than to F value. In fact, 
the fiux is maintained at high values {q > 0.9 x Qmax) for ^^,2/^^,1 ^ [0.9,1.1] and 
^2/^1 ^ [0.49, 2.53]. A 10% variation of Up leads to a comparable variation of the fiux 
value. The same fiux variation is obtained with a variation of F up to 150%. However, 
the asymmetry of g-behavior as a function of F is worth noting. In fact, the sign of 
F-variation affects strongly the variation of the fiux. Finally, the fiux sensitivity to ujp 
decreases for larger values of F. This is due to the fact that F controls the exchanged 
fiux spectral density peak width [10] : the larger F the larger the peak width which 
allows looser constraints on the peak position controlled by Up. 

4-2. Lorentz model 

First, we remind the dielectric permittivity expression according to this model t32j : 



ujq = ujto the transverse optical phonons circular frequency and F the damping factor. 



4-2.1. Identical media 

First, identical media are considered, medium 1 at 300 K and medium 2 at 299 K. 
The gap thickness between the two planes is 5 = 10 nm. Compared to Drude model, 
Lorentz model has an additional parameter, transverse optical phonons frequency cjq in 
this case. In this study, cjq = i^o^sic = 1.49 x 10-^^ (rad.s"-^) [32j is considered constant 
which reduces the problem to a two-parameter problem. Control parameters are ujlq 



e{uj) = Cc 
where uj^ uj\q — cJq 
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and r. Results will be presented as a function of and — . 
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Figure 8: Normalized NF RHF for 6oo = 1 (a-b) and eoo = 10 (c-d) obtained by exact 
(left column) and approximate (right column) calculations. See Table [s] for a summary 
of principal relevant results. 



Longitudinal phonons frequency and damping factor effect 

Figure [8] presents the normalized NF RHF exchanged by two semi-infinite planes which 
dielectric permittivities are modeled by Lorentz model for = 1 (8a 8b) and e^o = 10 



(8c 8d) obtained by exact (left column) and asymptotic calculations (right column). 
Principal relevant results of this figure, concerning the maximum position and value as 
well as calculation time, are summarized in Table |3| 

As for Drude model, we observe the existence of a maximum which is realized by a 
compromise between the phonons frequencies and the damping factor, i.e. between the 
peak position and width. For 6oo = 1 for instance, a maximal fiux density q^ax = 54529 
W.m"^ is observed at (^ = 1.42,^ = 1.9 x 10"^). Let us note that this (jmax 
value is almost five times lower than the value obtained with a Drude model at 6oo = 1 
{Qmax = 229336 W.m~^). In fact, Drude model is the Lorentz model limit when cjq 
goes zero. Thus, we observe, even though the detailed study of this parameter is not 
presented in the present paper, an increase of the maximal achievable fiux with a Lorentz 
model when ujq decreases. Furthermore, we observe that q^ax sensitivity to ujlq is much 
larger than to F. In fact, the fiux is kept at relatively high values {q > 0.9 x (jmax) 
with a relative variation of ujlo around ±12% {ujlo/^to ^ [1.25, 1.67]) versus a relative 
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variation of T up to +200% (T/loto e [7 x lO'^, 5.9 x lO'^]). 



Method 






UJTO 


r 


Qmax (W.m-2) 






t (s) 


E 


1 




1.42 


1.9 X 10"^ 


56896 


2000 


100 


2.8 X 10^ 


A 


1 




1.42 


1.9 X 10"^ 


56905 


10000 


100 


2.51 X 102 


E 


eoo,SiC = 


6.7 


1.24 


4.78 X 10-2 


14874 


4000 


100 


1.09 X 10^ 


A 


€oo,SiC = 


6.7 


1.24 


4.78 X 10-2 


14849 


10000 


100 


1.25 X 102 


E 


10 




1.22 


3.31 X 10-2 


10415 


3000 


100 


6.17 X 10^ 


A 


10 




1.22 


3.31 X 10-2 


10391 


10000 


100 


2.5 X 102 



Table 3: Optimal Lorentz model parameters for identical media plane-plane 
configuration. Values are obtained by exact (E) and approximate (A) calculations. 
^0 = ^^,SiC = 1-49 X 10-^^ rad.s~^ is kept constant, t is CPU time, n^^^ and n^^^r are the 
mesh points number for frequency and wave vector [uj and q) and control parameters 
(r/cjo and ujol/^o) discretization respectively. 



High frequency limit of the dielectric function effect 

As for Drude model and for the same reasons, we observe a decrease of q^ax when 6oo 
increases in addition to a shift of the maximum position to lower values of ujlo/^o 



Exact versus approximate calculation 

At this point, Lorentz model strongly contrasts with what was previously observed with 
Drude model giving very accurate asymptotic results for the maximal flux value as well 
as for its position. The position relative error is lower than 10~^. Similarly low relative 
error values are observed for the maximal flux value, except for the case e^o = 10 where 
it reaches 2.3 x 10~^. Indeed, the maximal flux relative error increases with eoo, i.e. 
when qjriax decreases. A part of this error is due to the omission of the propagative 
contribution in the asymptotic calculation. This contribution is almost constant for 
different values of e^o while the p-polarized evanescent contribution and the total flux 
decrease when increases. 

In spite of comparable accuracy, asymptotic calculations are still 1000 times faster than 
exact ones. In addition, the approximate method shows a better convergence. In fact, 
some numerical oscillations due to slow convergence can be observed on Figure 8c for 
small flux values. 

To understand the origin of the consistency of the two methods in the case of Lorentz 
model we proceed as done previously for Drude model and examine the transmission 
probability T^^ani^^Q) the integrand (q/qo) x ^evani^^Q) ^^^^ wave 

vector parallel component in equation ISl These two quantities are plotted in flgures 



9a and ^\ respectively. We first note that, compared to Drude model, modes are 



transmitted here in a much lower number which explains the lower fiux density values. 
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Figure 9: (a) Transmission probability of p-polarized evanescent modes T^yani^^Q) 
(b) the integrand (q/qo) x T^^ani^^ Q) the sum over the wave vector paraUel component 
in equation |8] for Lorentz model with e^o = 1- Lorentz model parameters are given in 
(Table |3) line 1). Only one case, the exact optimum in this case, is considered here since 
approximate calculations gave similar results with high accuracy. 



Second, transmitted modes mainly lay in the range q/q^ G [15, 150] if we consider 
^evani^^Q) — \' This Concentration of transmitted modes around medium and high q 
values is behind the high accuracy of the approximate method. Finally, the cutoff wave 
vector, considering the same criterion as in Drude model, is found for ^c/^g = 216, i.e. 
for qc ^ 1.8/5 which is of the order of 1/5. 

Case of silicon carbide at 300 K 

Finally, we consider the case of silicon carbide (SiC). This material has been extensively 
studied in NF radiative heat transfer literature for its strong surface phonon-polariton 
resonances around uj = 10^^ (rad.s~^). 




1.2 1.4 1.6 1.8 2 

Figure 10: Normalized NF RHF between two semi-infinite planes as a function of 
and ^ for e^^ = e^Q^siC = 6.7. SiC parameters values correspond to the point 
(1.24,6 X 10"^)[32J indicated by a cross. 



Figure 



10 



presents the normalized NF RHF exchanged by two semi-infinite planes 
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modeled by Lorentz model with e^o = ^oo.SiC = 6.7. With only 0.6 x (jmax^ SiC is far 
from approaching Lorentz model optimal performances unlike the case of HD-Si which 
parameters allowed flux values as high as 90% of the maximal RHF that can be obtained 
with a Drude model when = esi- Besides, this figure is obtained by calculations with 
higher resolution meshes, n^^q = 4000 in this case. Compared to figure 8c, this shows 
that numerical oscillations magnitude decreases slowly when the mesh points number 
n^^q of (cj, q) space increases. CPU time is however one order of magnitude larger than 



previously, i.e. than in figure 8c 



4.2.2. Non-identical media 

Now, consider two semi-infinite planes made of non identical materials. We will analyze 
two cases : (1) the case of SiC and a slightly different material (Figure 11a) (2) The 
case of the fictive material realizing the optimal performances with e^o = ^oo^SiC = 6.7 
(see Table [sj line 3) that we will note material 1 with a slightly different material 
(Figure lib). For both cases, ujq = ujto.SiC = 1.49 x 10-^^ rad.s~^ is constant. Control 



parameters are then (jJol/^ol.SiC and V /Vsic in the first case and (jJol/^ol,i and F/Fi 
in the second. 
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Figure 11: Normalized NF RHF exchanged by two semi-infinite planes of different non- 
identical materials which dielectric functions are modeled by Lorentz models. Two cases 
are considered : (a) SiC exchanging with another material, (b) the material maximizing 
the transfer between two identical media at = 6.7 (see Table [sj line 3) with another 
material. 



Unsurprisingly, the optimum is observed at (1, 1) in both cases, i.e. for identical media. 
We also observe a high sensitivity, more pronounced for SiC, of the fiux to ujlo- In fact, 
a relative variation of the order of 10~^ of ujlo around the point (1, 1) decreases the fiux 
below 0.9 X q^nax while a 8 x 10~^ relative variation of cj^o halves the fiux value, q is 
however much less sensitive to F since relative variations of this parameter in the range 
[—8%, 581%] maintains q > 0.9 x (jmax- The high asymmetry of this range around zero 
is due to the fact that F controls the imaginary part of the dielectric permittivity peak 
width and height. On the other hand, Im{e) controls both emission and absorption. 
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Thus, if r increases for the material exchanging with SiC, this material dielectric 
permittivity imaginary part peak will be wider and lower than SiC's. Therefore, all 
SiC modes will contribute to the transfer (the new peak is wider than SiC peak) , with a 
lower modes density though (the new peak is lower than SiC peak). Then again, when 
r decreases, the peak becomes narrower and higher than SiC's. All SiC modes do not 
contribute to the transfer anymore while contributing modes have the same density than 
in SiC-SiC system. 

Similar considerations can be made about he second case (Figure lib) with the only 
difference that material 1 damping factor is lower than SiC's. This implies a wider peak 
for Im{ei) which allows looser constraints on the peak position controlled by ujlq and 
cjq. For instance, a ±1.2% relative variation of ujlo keeps the flux higher than 0.9 x qmax- 
This value is two orders of magnitude higher than SiC's, even though it is still relatively 
small and restrictive in regard to the quality of materials that can be obtained with 
usual nano-materials deposition techniques. 



5. Conclusion 



In this work, a study of the effects of different parameters of usual materials local 
dielectric functions models (Drude and Lorentz) on NF RHF exchanged by two semi- 
infinite planes separated by a nanometric gap at room temperature is presented. For 
this purpose, exact and approximate (according to the asymptotic electrostatic limit 
approximation in the extreme near-field regime presented in ^5) calculations of the 
heat fiux were calculated. We then showed that the asymptotic approximation leads to 
highly accurate results, in particular for Lorentz model, with a calculation time at least 
one thousand times shorter than exact calculation time. Two particular materials usually 
considered for near-field heat transfer optimization were also considered : silicon carbide 
(SiC) and highly doped silicon (HD-Si). HD-Si reveals to be well adapted to this aim. 
In fact, it allows to reach 90% of maximal achievable heat fiux by a Drude model with 
= 11.8. It is however possible to overcome these performances by a metamaterial 
that would have a much lower value of the dielectric permittivity high frequency limit. 
On the other hand, Lorentz model in general, and SiC in particular, are not the best 
choice in order to maximize NF radiative heat transfer, at room temperature at least. 
In addition, SiC is particularly penalizing since its maximal performance is strongly 
dependent on the quality of used materials. Thus, very small discrepancies, of the 
order of 10~^, between the phonons frequencies of the two SiC samples would halve 
the maximal achievable radiative heat fiux. We also showed, for both models, that 
the maximal RHF is obtained when the two semi-infinite planes are made of identical 
materials. 

Finally, it is worth mentioning the mesoscopic description of NF radiative heat transfer 
recently developed p!Bl [T5] and which renews the understanding of this kind of 
transfer : radiative energy is transported through different modes which have different 
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transmission probabilities from one medium to the other. Total exchanged energy is 
then obtained by summing the energy of each mode weighed by the mode transmission 
probability. Maximizing the transfer reduces then to maximizing the transmission 
probability of the different modes. According to this idea but without performing a 
detailed optimization study, Ben-Abdallah and Joulain [16] derived with variations 
calculus a simple analytical condition on Fresnel reflexion coeflicients which allows, 
knowing medium 1, to determine Fresnel coeflicients of the second medium maximizing 
the transfer. It is then possible to determine the optical properties of both media. It 
would be interesting to implement this method and compare its results and performances 
to previously presented methods. 



6. Appendix : on the calculus of the cutoff wave-vector 

The cutoff wave vector Qc is the upper bound of q that would allow an accurate evaluation 
of the sum : 

revan(^,q) 

According to T^vani^^Q) (q/Qo) ^ ^evani^^Q) plc>ts (Figures [i] and [5] respectively. 
Drude model examples are considered to illustrate the method.), the cutoff wave vector 
Qc depends on the circular frequency u). At the present stage and for simplicity sake, 
we consider a constant cutoff wave vector Qc. According to the same flgures, the largest 
wave vectors participating to the transfer are observed for uj = uUp. Thus, the constant 
cutoff wave vector is to be determined at this frequency. Two families of criteria can be 
considered for Qc deflnition, whether this latter is based on the transmission coefficient 
^evani^^Q) ^hc trausmissiou coefficient weighted by the normalized wave vector, 



6.1. Transmission coefficient criterion 



A first criterion based on Tg^^^ (cj, q) can be considered. According to figure \12\ 
^evani^p^ ^) incrcascs with increasing ^ to reach its maximal value c::^ 1 at a certain wave 
vector qjriax and slowly goes to zero after that. Qc can be defined as the smallest wave 
vector larger than Qmax which separates transmitted modes from those with sufficiently 
small transmission probability, defined by an arbitrary threshold T^nin- Then Qc is defined 
by: 

I ^evanV^p-) qc) ~ ^min 

If we consider a threshold transmission probability Tj^in = 0.5 for example (this threshold 
value separates modes that are more likely to be transmitted from those who are not), 
this leads to qc = 2.73/5 and qc = 2.37/5 for the exact and the approximate optima 




respectively. A lower threshold of Tmin = 0.1 leads to Qc = 3.55/5 and Qc = 3.22 /S 
respectively. In all cases, Qc is of the order of 1/6. 



6.2. Weighted transmission coefficient criteria 



Two criteria have been considered : a direct one that can be directly verified on 
{q/Qo) X ^evani^^Q) color maps and an indirect one that needs further calculations. 



6.2.1. Full width at half maximum (FWHM) 
Figure 13 presents the integrand (q/qo) x r^^ani^p^Q) ^ 
If we consider that the monochromatic fiux density at = cj^ 



ujp. A peak is observed, 
is transmitted by modes 
with g in a wave vector range equal to the FWHM around the peak, then qc is the upper 
bound of the FWHM, i.e. the largest wave vector verifying : 



qc 

qo 



r{ujp,qc) 



T{ujp,q) 



(17) 



For Drude model for instance, this definition led to values of qc = 3/5 and qc = 2.7/5 
for the exact and approximate optima respectively. This criterion is interesting since it 
can be directly verified on ^^{uj^ q) color maps and was used to obtain results reported 
in the present paper. 



6.2.2. Fractional monochromatic flux 

A second criterion, more complicated to implement though since it can not be read on 
-^T{uj^q) maps and needs an integral calculation, considers the monochromatic radiative 
heat fiux fraction transmitted in a certain wave vector range. 




We can then define Qc as the wave vector that verifies (see Figure 14) : 
q q 

-^r{ujp,q)dq = Xf X I -^r{u)p,q)dq (18) 
c- % % 
where Xi^^ is the monochromatic fiux density transmitted fraction. This second criterion 
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Figure 14: Cutoff wave vector definition according to the fractional monochromatic NF 
RHF density carried by modes in the range g G [go, Qc\- 



is expected to be more accurate than previously presented ones since it provides a 
rigorous quantitative information, xf in this case. For xf = 0.95 for example, it leads 
to qc = 3.6/5 and qc — 3/S for the exact and the approximate optima respectively. 
The correction compared to the previous criterion results varies from 20% to 10% 
respectively. However, qc is still of the order of 1/5. 
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